Astronomy &: Astrophysics manuscript no. 2416 


February 2, 2008 


(DOI: will be inserted by hand later) 





The CMB temperature power spectrum from an improved 

analysis of the Archeops data 

M. Tristram^ G. Patanchon^, J. F. Maci'as-Perez\ P. Ade^, A. Amblard"', R. Ansari^, E. Aubourg^''^, A. 
Benoit^, J.-Ph. Bernard^, A. Blanchard^°, J. J. Bock^^-^^, F. R. Bouchet^^, A. Bourrachot^, P. Camus^, 
J.-F. Cardoso^^, F. Couchot^ P. de Bernardis^^ J. Delabrouille^, F.-X. Desert^^^ M. Douspis^°, L. 
]Q ; Dumoulini^ Ph. Filliatrei^'^ P. Fosalba^^ M. Giard^, Y. Giraud-Heraud^ R. Gispert^Of*, L. Guglielmi^ 

J.-Ch. Hamilton^^, S. Hanany^^, S. Henrot-Versille^ , J. Kaplan^, G. Lagache^°, J.-M. Lamarre^^, A. E. 
^ ; Langei\ K. Madet^ B. MafFei^, Ch. Magneville*^' S. Masil^ F. Mayet\ F. Nati^^ O. Perdereau^ S. 
. '. Plaszczynski^, M. Piat^, N. Ponthieu^^, S. Prunet^^, C. Renault\ C. Rosset"^, D. Santos\ D. Vibert^^, and 
d' D. Yvon^ 



o 



X 



Laboratoire de Physique Subatomique et de Cosmologie, 53 Avenue des Martyrs, 38026 Grenoble Cedex, France 
Department of Physics & Astronomy, University of British Columbia, Vancouver, Canada 
pi ; ^ Cardiff University, Physics Department, PO Box 913, 5, The Parade, Cardiff, CF24 3YB, UK 

^ '. * University of Cahfornia, Berkeley, Dept. of Astronomy, 601 Campbell HaU, Berkeley, CA 94720-3411, U.S.A. 

(yT) ' ^ Laboratoire de I'Accelerateur Lineaire, BP 34, Campus Orsay, 91898 Orsay Cedex, France 

fT^ [ ® CEA-CE Saclay, DAPNIA, Service de Physique des Particules, Bat 141, F-91191 Gif sur Yvette Cedex, France 

\0 • APC, 11 pi. M. Berthelot, F-75231 Paris Cedex 5, France 

Centre de Recherche sur les Tres Basses Temperatures, BP166, 38042 Grenoble Cedex 9, France 
Centre d'Etude Spatiale des Rayonnements, BP 4346, 31028 Toulouse Cedex 4, France 
nI ■ ^'^ Laboratoire d'Astrophysique de Tarbes Toulouse, 14 Avenue E. Belin, 31400 Toulouse, France 

California Institute of Technology, 105-24 Caltech, 1201 East California Blvd, Pasadena CA 91125, USA 
Jet Propulsion Laboratory, 4800 Oak Grove Drive, Pasadena, California 91109, USA 
Institut d'Astrophysique de Paris, 98bis, Boulevard Arago, 75014 Paris, France 
JL I " CNRS-ENST 46, rue Barrauh, 75634 Paris, France 

< I ■ Gruppo di Cosmologia Sperimentale, Dipart. di Fisica, Univ. 'La Sapienza', P. A. Moro, 2, 00185 Roma, Italy 

' Laboratoire d'Astrophysique, Obs. de Grenoble, BP 53, 38041 Grenoble Cedex 9, France 



11 



5j . " CSNSM-IN2P3, Bt 108, 91405 Orsay Campus, France 

CEA-CE Saclay, DAPNIA, Service d'Astrophysique, Bat 709, F-91191 Gif sur Yvette Cedex, France 



19 
20 



Institute for Astronomy, University of Hawaii, 2680 Woodlawn Dr, Honolulu, HI 96822, USA 
Institut d'Astrophysique Spatiale, Bat. 121, Universite Paris XI, 91405 Orsay Cedex, France 



5_j . 21 LPNHE, Universites Paris VI et Paris VII, 4 place Jussieu, Tour 33, 75252 Paris Cedex 05, France 

I School of Physics and Astronomy, 116 Church St. S.E., University of Minnesota, Minneapolis MN 55455, USA 

LERMA, Observatoire de Paris, 61 Av. de I'Observatoire, 75014 Paris, France 

February 16, 2005 

Abstract. We present improved results on the measurement of the angular power spectrum of the Cosmic 
Microwave Background (CMB) temperature anisotropics using the data from the last Archeops flight. This 
refined analysis is obtained by using the 6 most sensitive photometric pixels in the CMB bands centered at 143 
and 217 GHz and 20 % of the sky, mostly clear of foregrounds. Using two different cross-correlation methods, 
we obtain very similar results for the angular power spectrum. Consistency checks are performed to test the 
robustness of these results paying particular attention to the foreground contamination level which remains well 
below the statistical uncertainties. The multipole range from ^ = 10 to ^ = 700 is covered with 25 bins, confirming 
strong evidence for a plateau at large angular scales (the Sachs-Wolfe plateau) followed by two acoustic peaks 
centered around £ = 220 and £ = 550 respectively. These data provide an independent confirmation, obtained at 
different frequencies, of the WMAP first year results. 
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1. Introduction 

Observations of the Cosmic Microwave Background 
(CMB) temperature anisotropics provide answers to 
fundamental questions in cosmology. The experimental 
determination of the CMB temperature angular power 
spectrum (|Netterfield et al.(1997)] ^lillcr et al. 1999, 

IHanany et al. 200d| 
INette rfield et al. 20021 
[Sievcrs et al. 200ji 
[B^it et a I. 2nn8al 
iB^ats et al. 20041 



IdeBernardis et al. 2000r 
ILce et a l. 200T1 
[Halvcrson et al. 20021 
IRuhino-M artin et aJ. 20081 
[Hinshaw_ 



iRcadh cad et al. 20041 ILeitch et al. 2004|l leads to im- 
portant insights into the composition and evolution of 
the Universe. Most notable are the conclusions that 
the geometry of space is essentially flat, the measure- 
ments are consistent with the inflationary paradigm and 
the Universe is dominated by unknown forms of dark 
energy and dark matter 
IMacfas-Perez et al. 2000l 



( Lineweaver et al.(1997) 
IBenoit et al. 2003b 



IDouspis et al. 2003|[^ergel et al. 2003| ). 

Archeops^ was designed to obtain a large sky cov- 
erage of CMB temperature anisotropics in a single bal- 
loon flight at millimeter and submillimeter wavelengths. 
Archeops is a precursor to the Planck HFI instrument 
IjLamarre et al. 207)3|l . using the same optical design and 



the same technology for the detectors, spider-web bolome- 
ters, and their cooling, 0.1 K dilution fridge. The instru- 
ment consists of a 1.5 m aperture diameter telescope and 
an array of 21 photometric pixels operating at 4 frequency 
bands centered at 143, 217, 353 and 545 GHz. The data 
were taken during the Arctic night of February 7, 2002 
after the instrument was launched by CNES from the 
Esrange base near Kiruna (Sweden). The entire data set 
covers ^ 30% of the sky. 

The Archeops initial analysis IjBenoit et al. 2003a|l 
- hereafter Paper I - presented for the first time mea- 
surements from large angular scales to beyond the first 
acoustic peak {£ — 15 — 350). A few months later, the first 
year WMAP results IjBennett et al. 2003f) confirmed the 
previous measurements and significantly reduced the error 
bars on scales down to the second acoustic peak. 

This paper presents a second and more refined analy- 
sis of the Archeops data. With respect to Paper I, ma- 
jor improvements on the timeline processing, the map- 
making, the beam modeling and the foreground removal 
were achieved. Further, new power spectrum estimation 
methods based mainly on the cross power spectra between 
different detectors maps are used to reduce the contri- 
bution from correlated noise and systematic effects. This 
essentially allows us to increase the number of detectors 
considered (from two for Paper I to six for this analysis) 
and to cover a larger fraction of clean sky (12 % in Paper 
I, 20 % in this paper). These developments lead to a bet- 
ter sampling and a larger range in multipole space with 
an improved accuracy. 



The paper is organized as follows: Sect. |21 summa- 
rizes the processing on the TOIs (Time Ordered Data) 
with an emphasis on changes and improvements with re- 
spect to Paper I. Section 13 describes the methods, Xspect 
and SMICA, used for the estimation of the CMB angular 
power spectrum from observed emission maps. The esti- 
mation of the Archeops CMB angular power spectrum is 
presented in Sect.^ Consistency checks on the data and 
the contribution from systematics to the Archeops CMB 
angular power spectrum are discussed in Sect. El A sim- 
ple comparison with the best-fit cosmological model pro- 
vided by the WMAP team ( |Spergel etal. 2003| ) is shown 
at the power spectrum level. However, we postpone to a 
forthcoming paper the comparison of this dataset to the 
WMAP data and other datasets at the map level. 

2. Observations and Data processing 

The Archeops experiment is described in details in com- 
panion papers. Instrument and data processing are de- 
tailed in lMacias-Perez et al. 2005|l while the in-flight per- 
formances are summarized in lMadet et al. 2004a] In the 
following subsections, only key points on the data process- 
ing are summarized and we then focus on refinements im- 
plemented for the present analysis, as compared to Paper 



2.1. Observations and standard data processing 

The instrument contains a bolomctric array of 21 pho- 
tometric pixels, each one being made of cold optics con- 
sisting of an assembly of back-to-back horns, filters and 
lenses, and of a 100 mK bolometer, which operate at fre- 
quency bands centered at 143 GHz (8 pixels), 217 GHz 
(6), 353 GHz (6=3 polarized pairs) and 545 GHz (1). The 
two low frequencies are dedicated to CMB studies while 
high frequency bands are sensitive essentially to interstel- 
lar dust and atmospheric emission. The focal plane is made 
of 21 spider-web bolometers and some thermometers and 
is maintained at a temperature of ~95 mK by a ^He- 
^He open-circuit dilution cryostat. Observations are car- 
ried out by spinning the payload around its vertical axis 
at 2 rpm. Thus the telescope produces circular scans at a 
fixed elevation of ^ 41 deg. Observations of a single night 
cover a large fraction of the sky as the circular scans drift 
across the sky due to the rotation of the Earth and the 
gondola trajectory. 

The Archeops experiment was launched on February 
7, 2002 by the CNES^ from the balloon base in Esrange, 
near Kiruna, Sweden, 68°N, 20°E. The night-time scien- 
tific observations span 11 hours of integration. The point- 
ing reconstruction, with rms error better than 1 arcmin., 
is performed using data from a bore-sight mounted opti- 
cal star sensor. Each photometric pixel offset is deduced 
from Jupiter observations. 



^ see jhttp : //www . archeops . org^ 



^ Centre National d'Etudes Spatiales, the French national 
space agency 
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Corrupted data (including glitches) in the Time 
Ordered Information (TOI), representing less than 1.5% of 
the full data set, are flagged. Low frequency drifts on the 
data generally correlated to house-keeping data (altitude, 
attitude, temperatures, CMB dipole) are removed using 
the latter as templates. Furthermore, a high frequency 
decorrelation is performed in few chosen time frequency 
intervals of ~1 Hz width to remove some bursts of non- 
stationary high-frequency noise localised in time and in 
frequency. The corrected timelines are then deconvolved 
from the bolometer time constant and the flagged cor- 
rupted data are replaced by a realization of noise (which 
is not projected onto the maps in the map-making step). 
Finally, low time frequency atmospheric residuals are sub- 
tracted using a destriping procedure which slightly filters 
out the sky signal to a maximum of 5% (see the red curve 
on Fig 01). This effect is corrected for when computing the 
CMB angular power spectrum as discussed in Sect. 13.21 

The CMB dipole is the prime calibrator of the instru- 
ment. The absolute calibration error against the dipole as 
measured by COBE/DMR IjFixsen et al. 1994|l and con- 
firmed by WMAP (jBennett et al. 2003|l is estimated to 
be 4% and 8% in temperature at 143 GHz and 217 GHz 
respectively. These errors are dominated by systematic ef- 
fects. 

As Jupiter is a point-source at the Archeops resolu- 
tion, local maps of Jupiter allow us to estimate the time 
constant of the bolometers and the main beam shape. This 
is performed using the two Jupiter observation windows. 
While the 143 GHz detector beams are mostly elliptical, 
the 217 GHz ones are rather irregular (multi-mode horns). 
The typical FWHM of the beams is about 12 arcmin. Two 
Saturn crossings allowed cross-checks on the time con- 
stants and beams. 

2.2. Removal of Galactic and atmospheric foreground 
emissions 

The Archeops cleaned TOIs at 143 and 217 GHz are 
contaminated by atmospheric residuals coming mostly 
from the inhomogeneous ozone emission. This contributes 
mainly at frequencies lower than 2 Hz in the timeline and 
follows approximatively a law in antenna temperature. 
Therefore atmospheric emission is much more important 
at the high Archeops frequencies (353 and 545 GHz). In 
the same way, at the Archeops CMB frequencies (143 
and 217 GHz) the Galactic dust emission also contami- 
nates the estimation of the CMB angular power spectrum 
even at intermediate Galactic latitudes. Dust emission, 
which presents a modified black-body spectrum at about 
17 K with an emissivity of about z^^, dominates the CMB 
at high frequencies and therefore the 353 and 545 GHz 
channels can be used to monitor it. To suppress both resid- 
ual dust and atmospheric signals, the data are decorre- 
lated using a linear combination of the high frequency pho- 
tometric pixels (353 and 545 GHz) and of synthetic dust 
timelines. These are constructed from the extrapolation 




Fig. 1. Galactic mask (dark blue lane) applied to the 
Archeops coverage (annular green region). The CMB 
mask is obtained by requesting the SFD brightness at 
353 GHz to be < 0.5 MJy.sr'^ The MoUweide projection 
of the celestial sphere is in Galactic coordinates centered 
on the Galactic anti-center. Gridding on the full sky map 
is by 30 degree steps. The CMB analysis includes 20 % of 
the sky (dark green area) while Archeops covers ~30 % 
of the sky. The previous analysis only covered 12 % of the 
sky above the 30 degree Northern parallel. 

of IRAS and COBE observations in the far infrared do- 
main ( ISchlegel et"al. 1998, Finkbciner et al. 19 99) to the 
Archeops frequencies. We actually construct a synthetic 
dust template for the considered CMB bolometers and 
also for the high frequency bolometers so that we can take 
into account simultaneously in such a model both types 
of frequency behaviors. 

As the decorrelation is not perfect in the Galactic 
plane, a Galactic mask is then applied to the Archeops 
maps for determination of the CMB power spectrum. This 
mask is deduced from a Galactic dust emission model 
l lSchlegel et al. 1998| IFinkbeiner et al. 1999|l at 353 GHz. 
The Galactic plane and the Taurus region are efficiently 
masked by considering only regions with a brightness 
< 0.5 MJy.sr"^. Applying this mask, the CMB maps de- 
rived from the Archeops data cover 20 % of the sky 
sampled by ~ 100,000 pixels of 7 arcmin. (HEALPix 
nside — 512). Figure ^ presents the Archeops coverage 
to which we have superimposed the Galactic mask. Only 
the Northern part above 30 degree was used in Paper I. 

2.3. Map-making 

The noise power spectrum of the Archeops TOIs is 
nearly flat with increasing power at very low time fre- 
quencies due to residuals from atmospheric noise, and at 
very high time frequencies due to the deconvolution from 
the bolometer time constants. To cope with these two fea- 
tures on the Archeops noise we have used an optimal 
(i.e. it achieves least square error on pixelised map) pro- 
cedure called MIRAGE \Yvon fc Mayet 2004| ) to produce 
maps for each of the detectors. 

MIRAGE is based on a two-phase iterative algorithm, 
involving optimal map-making together with low fre- 
quency drift removal and Butterworth high-pass filtering. 
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Fig. 2. Archeops map of the CMB sky in Galactic coor- 
dinates centered on the Galactic anti-center after smooth- 
ing with a 30 arcmin. Gaussian. A patch of the sky of 
30x30 deg., with high redundancy and centered on {I, b) 
= (195,45) degrees is zoomed up. Gridding on the full sky 
map is by 30 degree steps, gridding on the zoomed patch 
is 5 deg. The Galaxy is masked as described above. 

A conjugate gradient method is used for resolving the lin- 
ear system. A very convenient feature of MIRAGE is that 
it handles classic experimental issues, such as corrupted 
samples in the data stream, bright sources and Galaxy 
ringing effects in the filtering and in the calculation of the 
noise correlation matrix. 

Maps are computed with 7 arcmin. pixels (HEALPix 
nside — 512) for each absolutely calibrated detector with 
their data time band-passed between 0.1 and 38 Hz. This 
corresponds to about 90 deg. and 20 arcmin. scales, re- 
spectively. The high-pass filter removes remaining atmo- 
spheric and Galactic contamination. The low-pass filter 
suppresses non-stationary high frequency noise. 

About two thirds of the Archeops sky are observed 
with ~20 to 60 samples per bolometer and per square 
degree and one third with a higher redundancy, about 75 
samples per bolometer and per square degree. For illustra- 
tion. Fig. 121 shows a map obtained from a weighted linear 
combination of the maps of each of the six most sensitive 
Archeops detectors. This map is smoothed with a 30 ar- 
cmin. gaussian beam and has a typical rms noise of 50 /iK 
per 30 arcmin pixel. 

3. Power spectrum estimation 

In this section, we present three methods, Xspect 
UTristram et al. 2004|l . SMICA UPatanchon 20038 
and power spectrum on the rings (Fm hereafter) 
IjAnsari et al. 2003)) used for the determination of the 
angular power spectrum of the CMB temperature 
anisotropics with the Archeops data. Beforehand we 
detail the procedure we use to correct from beam smooth- 
ing and filtering effects as well as from inhomogeneous 
coverage. 

We have thoroughly probed Xspect and SMICA with 
simulations which are described below. Results from both 
methods are included in this paper to cross validate the 



final results. The F^ method is provided here to illus- 
trate its potential in the estimation of the angular power 
spectrum directly from ring data and is more suitable to 
Planck-like data. 

Xspect and SMICA are based on the so-called 'pseudo- 
Ce's estimators (|Peebles 19731 |Szapudi et al. 2001[ 
IHivon et al. 2002|l which directly compute the pseudo 
power spectrum from the spherical harmonics decompo- 
sition of the maps. These spectra are then corrected from 
the sky coverage, beam smoothing, data filtering, pixel 
weighting and noise biases. 

A pseudo power spectrum Di is linked to the true 
power spectrum C'i by 



D, ^Y,Mu,plBlTi,{Ci>) -f {Ni) 



(1) 



where Mui is the mode-mode coupling matrix, Bi is the 
beam transfer function describing the beam smoothing ef- 
fect. Pi is the transfer function of the pixelization scheme 
of the map describing the effect of smoothing due to the 
finite pixel size and geometry, Ti is an effective transfer 
function that represents any filtering applied to the time 
ordered data, and {Ni) is the noise power spectrum. 

In the following, the Mui matrix describes the mode- 
mode coupling resulting from the incomplete sky cover- 
age and the weighting applied to the sky maps. We take 
into account the pi pixel transfer function due to the 
smoothing effect induced by the finite size of the map 
pixels. This function is provided in the HEALPix package 
l|(;orski et al. 199911. 

3.1. Beam smoothing effect 

Most of the beams of the Archeops detectors have been 
measured on Jupiter to be elliptical. A few of them are 
irregular. Therefore, the effective beam transfer function 
must be carefully estimated for each bolometer. The beam 
transfer functions are computed from simulations using 
the Asymfast method detailed in IjTristram et al. 2004|l . 
This method is based on the decomposition of the beam 
into a sum of Gaussians for which convolution is easy 
in the spherical harmonic space (up to 12 Gaussians are 
used here). This allows us to deal with asymmetric beam 
patterns using the scanning strategy of the instrument. 
Figure shows the beam transfer function for each of 
the Archeops detectors used in this analysis. They are 
estimated with a Monte-Carlo of 100 Asymfast simula- 
tions per bolometer. The beam transfer functions for the 
143 GHz detectors are very similar and close to circular 
Gaussian. The 217 GHz detector beams are larger and 
more irregular, and smear-out more the high multipoles. 

The Asymfast method produces negligible (< 0.1 %) 
statistical uncertainties on the Bi estimation. However, as 
the beam patterns have been measured on Jupiter maps 
they may differ from the effective beams on the CMB 
anisotropics. This comes mainly from uncertainties on the 
electromagnetic spectral dependence, far-side lobes, base- 
line subtraction and time constants, each of which esti- 
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Fig. 3. Beam transfer functions of the six most sensitive 
Archeops detectors computed using the Asymfast beam 
description. 



Fig. 4. Archeops transfer functions: Fi filter function 
(in blue), Di destriping transfer function (in red) and Ti 
total Archeops transfer function (in black). 



mated to be lower than 5%. For such systematics it is dif- 
ficult to estimate their impact on the beam transfer func- 
tion. As an illustration, we give conservative upper limits 
on the Be uncertainties by taking, as 1-sigma level error, 
a third of the difference between resulting transfer func- 
tion from elliptical beams IjFosalba et al. 2002|l and that 
from the Asymfast decomposition in multiple Gaussians. 
Fig.|n|shows the uncertainties on the CgS due to the beam 
transfer function uncertainties. They are well below the 
statistical error bars. 

3.2. Filtering and inhomogeneous coverage effects 

Filtering leads to a preferred direction on the sky (the 
scanning direction) and so the assumption of isotropic 
temperature fluctuations implicitly done in Eq. ^ is not 
valid any more. However, to a first approximation, the 
bias on the CMB power spectrum due to the filtering of 
the time ordered data can be accounted for in the spherical 
harmonic space through the Ti transfer function. 

For this analysis we have performed two types of fil- 
tering associated with the destriping of the data discussed 
in Sect. 12.11 and with the band-pass filter applied to the 
data on the map making procedure. 

The band-pass filter function Fg is computed from 
100 simulations of the CMB sky. The simulated maps 
are converted into timelines using Archeops pointing. 
These timelines are then filtered as the Archeops data. 
Subsequently, they are projected onto maps and the power 
spectrum of those is compared to the power spectrum ob- 
tained from maps of the same but unfiltered timelines. 



Figure 01 shows in blue the band-pass filter function. 
It reaches 65 % at € = 10 and remains above 85 % in the 
multipole range [25-700]. In our analysis, all bolometers 
are identically filtered and the difference between their 
pointing vectors is very small as these bolometers are dis- 
tributed onto two rows separated by only ^30 arcmin in 
the focal plane. We therefore assume an identical Fg func- 
tion for all detectors. Uncertainties on the estimation of 
the Fe function are derived from the dispersion of the sim- 
ulations. 

The transfer function associated with the destriping, 
De, has been computed using simulations and is shown in 
red on Fig. 01 The accurate determination of this function 
is difficult because the destriping procedure is non linear 
and CPU intensive. Thus, in order to be very conservative, 
we choose to take a third of the estimate of the function 
itself as the systematic error for it. 

The total transfer function used for the Archeops 
pipeline Ti = Fg x Dg is plotted in black on Fig. 0] The 
uncertainties on the final power spectrum due to the errors 
on the Tg function are represented on Fig. 

3.3. Xspect 

The Archeops angular power spectrum has been com- 
puted using an extension of the 'pseudo-C^' method to 
cross power spectra called Xspect IjTristram et al. 2004|) . 
Assuming no noise cross-correlation between different de- 
tectors, the noise term in Eq. vanishes and each cross 
power spectrum, A ^ B, is an unbiased estimate of the 
Cgs. Pseudo cross power spectra can be easily corrected 
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from inhomogeneous sky coverage, beam smoothing and 
filtering effects by extending Eq. n]into: 

5t^^J2 Mf,?plB^Bf,T,{Cr) (2) 

I' 

wliere tlie beam transfer functions Bi for each bolome- 
ter and the transfert function Tf are those previously de- 
scribed. The mode-mode coupling kernel M^^P is com- 
puted for each cross power spectra from the cross-power 
spectrum of the weighted masks. For the noise weight- 
ing scheme we consider a different noise weighted mask 
for each Archeops detector. This mask is constructed by 
multiplying the mask on Fig. ^ by the inverse of the noise 
variance on each pixel and is convolved by a 30 arcmin. 
Gaussian. 

After correction, all cross power spectra (C^^") are 

combined into a single estimate of the power spectrum, C'e, 
by weighted averaging assuming the correlation between 
multipoles to be negligible. This last assumption is not 
completely true, as we can see some correlation at low 
multipoles on Fig. [7\ Thus the estimate is not completely 
optimal but no measurable bias has been found in tests 
of Xspect on realistic simulations of Archeops data sets. 
Analytical estimates of the covariance matrix and of the 
error bars in the power spectrum are also given. 

Xspect is designed to estimate both the angular power 
spectrum and its error bars even with incomplete sky cov- 
erage and mask inhomogeneities, as is the case with the 
present Archeops data. The approach has been validated 
with simulations including realistic noise and CMB tem- 
perature anisotropics. The noise timelines are simulated 
from an estimation of the Fourier power spectrum of the 
noise IjAmblard &: Hamilton 2004) for each of the pho- 
tometric pixels. The CMB signal is simulated using the 
HEALPix software from the Archeops best-fit ACDM 
model (^Benoit et al. 2003b|l convolved by the beam trans- 
fer function. Signal and noise are added into a single time- 
line which is filtered as the Archeops data and projected 
on the sky using the Archeops pointing. 

Three sets of 1000 simulations have been computed for 
sky maps with HEALPix resolution nside — 512 : a first 
one using an uniform weighting, a second one using a noise 
weighting scheme, and a third one with no noise added. 
Simulations were performed using the same optimal map- 
making method fYvon fc Mayet 2004| ) as the one used for 
the data. 

From these simulations we have found that there is no 
bias at the 1% level in the estimation of the power spec- 
trum. The analytical error bars provided by Xspect are 
also found to be above the standard deviation in the sim- 
ulations by less than 10% and with a rms of 7%. Moreover, 
the noise contribution to the error bars on the simulated 
data and the Archeops data are in agreement within 
5%. Hereafter, we will use the analytical estimates pro- 
vided by Xspect for the error bars of the Archeops an- 
gular power spectrum excluding the sample variance con- 
tribution. The latter is computed from the dispersion of 



the simulations without noise and is added up to obtain 
the final error bars on the CMB angular power spectrum. 
Therefore, the sample variance contribution to the error 
bars is given by the best-fit Archeops model described 
inlRenoit et al. 2008hl 

As mentioned ealier, an improvement of about 10% on 
the error bars is obtained by using uniform weighting at 
low multipoles and a noise weighting scheme at high mul- 
tipoles. Thus, in the following all power spectra presented 
are computed using uniform weighting up to ^ < 260 and 
using a noise weighting scheme for £ > 260. 

3.4. SMICA 

Using the filtering and beam transfer functions as well 
as the masks described in Sect. 13.11 and 13.21 we pro- 
cess the Archeops maps with a different estimation 
method of the CMB angular power spectrum: SMICA 
(Spectral Matching Independent Component Analysis) 
UPatanchon 2003|l . 

A specificity of SMICA is its ability to estimate jointly 
the power spectra of several underlying components (in- 
cluding noise) assuming that the observed sky is a linear 
combination of components. In spherical harmonic space 
and in a matrix form, the model is : 

X£m = A Se,n + n^m (3) 

where xim is a vector of spherical harmonics coefficients 
of the observed maps for each of the considered detec- 
tors; A is the Nd (number of detectors) x Nc (number 
of components) mixing matrix which defines the ampli- 
tude of the different components in each observed map. 
The coefficients of A are related to the electromagnetic 
spectra of the components and to the relative calibration 
between detectors. The spherical harmonic coefficients of 
the components and noise are stored in vectors sim and 
ngm. 

SMICA is based on matching empirical auto- and cross 
spectra to their expected forms, as predicted by model (jj)) 
and by the statistical assumption of decorrelation between 
components. The mismatch is measured by a measure 
of divergence between the measured and modeled spec- 
tra which stems from the likelihood of a Gaussian sta- 
tionary model. The adjustable parameters are: the power 
spectrum of each of the components (including CMB and 
noise) as well as the mixing matrix A. A complete de- 
scription of SMICA is given in IDelabroiiille et al. 20031 
l( Cardoso et a,l."2002l IPata,nchori 20031 

In the specific case of Archeops, spectral statistics are 
formed as follows. The spherical harmonic coefficients xim 
are computed on the sky region which is common to all de- 
tectors using two different weighting schemes. For i < 260, 
pixels are uniformly weighted. For £ > 260, pixels are 
weighted proportionally to the number of data samples per 
pixel for the best detector. Band-averaged pseudo auto- 
and cross-power spectra are formed from these Xim and 
corrected for beam smoothing. If Q bands are used, we 
obtain in this manner a set of Q spectral matrices Rq 
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{q — 1, . . . , Q), each of size Nd x Nd- Next, we choose 
which parameters should be be estimated (power spectra 
for CMB and possibly other components, all or parts of 
the coefficients, noise levels), collect all these parameters 

into a vector and denote Rq{9) = (^^g^ the expected 
value of the spectral matrices for a given value of 9 (this 
is easily computed from model (O). The SMICA algo- 
rithm estimates the unknown parameters by minimizing 
the spectral mismatch 

(/.(0) = ^u;,X(i?„i?,(0)) (4) 
q 

where Wq is the number of independent aim iu the 
qth spectral band and where the mismatch measure 
K{-,-) between two positive matrices is defined as 
K{Ma,Mb) = \ (trace(Af,M,-i) - logdetM.M"! -N^) 

(with this choice, the estimated parameter 9 — 
argmin(/>(0) is a maximum likelihood estimate as shown 
in IDelabrouille et al. 2003|l . The resulting estimated 
power spectra are then corrected from partial coverage and 
filtering effects using the MASTER formalism described in 
Sect.lO 

In order to evaluate error bars and possible biases, 
we have performed 500 realistic simulations of Archeops 
data. The data model includes synthetic CMB emission 
(observed with the same scanning strategy as used by 
Archeops) and noise for each detector. Application of 
SMICA to these simulated data has not shown any mea- 
surable bias. 

Error bars for the estimated power spectra can also be 
obtained analytically from the Fisher information matrix. 
They have been compared to the dispersion found in the 
Monte-Carlo simulations. Analytic error bars on the CMB 
power spectrum are found to be slightly underestimated 
(about 10% on average). In the following, we use the ana- 
lytic error bars corrected from the factor measured in the 
simulations. 

3.5. CMB power spectrum on the rings 

A third approach based on one-dimensional prop- 
erties of the CMB inhomogeneities on rings has 
been performed on Archeops data IjAnsari et al. 20031 
|PIaszczynski fc: Couchot 2003) . It has been made possi- 
ble by the Archeops sky scanning strategy, which scans 
quasi circles on the sky. The fact that we directly use TGI 
information with no requirement of projection on maps 
of the sky makes this method complementary to the two 
previous ones. 

Tm is defined as the Fourier power spectrum of the 
signal on a sky ring. For a ring of colatitude Q, the rela- 
tion between rm(0) and the Ce l|Delabrouille et al. 199811 
follows: 

oc 

r™(e) = cmBj{v)lm{Q) . (5) 

t=\m\ 



where is the transfer function for the destriping and 
filtering, Bg, is the beam transfer function and ('P)^^ are 
the Legendre polynomials. 

Rings are built for each bolometer from the TOIs by 
using the pointing information. They are then analysed by 
pairs. For each ring pair (i — l,i), whenever measurements 
taken at the same angular phase are separated on the 
sky by less than 0.1 degree, we define a "signal" Si{(j)) 
and a "noise" Ni{<j)) as respectively the half sum and half 
difference of the measurements from each rings. 

Once these quantities are computed ring per ring, we 
analyse S and N in two ways. On the first hand, we com- 
pute the difference of the mean values of their Fourier 
spectra (that we call the analysis). On the other hand, 
the average of the auto-correlation functions for each pair 
is computed and then Fourier transformed to obtain the 
Tm power spectrum. In both cases a Galactic mask similar 
to that described in Sect. 12. 21 is apphed. In addition since 
the autocorrelation approach needs all the low frequency 
drifts to be properly removed, we apply a cross-scan de- 
striping HBourrachot 2004|l . Since the noise directly pops- 
up from the data themselves, no simulation is needed in 
these approaches. 

The error bars on the r,„ power spectrum are com- 
puted from the dispersion on the Fourier transform across 
rings and then propagated to obtain the uncertainties on 
the angular power spectrum. 

4. Main results 

The analysis presented in this paper uses the six most sen- 
sitive Archeops bolometers, four at 143 GHz and two at 
217 GHz with instantaneous sensitivities ranging from 93 
to 207 /iKcA/B S2 . Note that those instantaneous sensitivi- 
ties are better, by a factor of at least five, than those of the 
WMAP satellite mission detectors (Bennet t et al. 20031 
and a factor 2 to 4 worse than the nominal ones expected 
for the Planck-HFI instrument. We consider 20 % of 
the sky by applying the Galactic mask presented in Sect. 

Table^presents the angular power spectrum measured 
by Archeops. Results for the Xspect and SMICA meth- 
ods are both given as they are based on different assump- 
tions on the data model. 

4.1. Archeops temperature angular power spectrum 
using Xspect 

Figure [5l shows the Archeops CMB angular power spec- 
trum obtained using the Xspect method for two inter- 
twined binnings (blue and red) . These binnings correspond 
to two sets of overlapping and shifted window functions 
which lead to two non-independent estimates of the CMB 
angular power spectrum. A mix of logarithmic and lin- 
ear scales in multipole space is presented to improve the 
readibility of the figure both on the Sachs-Wolfe plateau 
and on the first two acoustic-peaks clearly detected by 
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Fig. 5. Archeops temperature angular power spectrum 
obtained using the Xspect method. A mixing of log-hnear 
scales is presented to improve the readibility of the fig- 
ure both on the Sachs-Wolfe plateau and on the acoustic 
peaks regions. Two intertwined and therefore not indepen- 
dent binnings (red and blue) are represented. 

Archeops. Two different weighting schemes are com- 
bined to produce the smallest error bars. At low (. mul- 
tipoles a uniform weighting is preferred whereas for high 
£s the sky maps for each detector are noise weighted by 
using Wp^d = ^/o'p d "where ^ is the variance of the pixel 
p of the sky map from the detector d. The two schemes 
yield identical results around the mixing point, £ ~ 260 
and they are joined in order to minimize the final error 
bars. 

Figure shows a detailed description of the statisti- 
cal error bars (in black) on the Archeops angular power 
spectrum in terms of sample variance (in cyan) and instru- 
mental noise (in red). Sample variance is deduced from the 
set of simulation without noise. It corresponds to the un- 
certainty on the model that is induced by the fact that 
we can only look at a part of one realisation of the sky. 
Sample variance dominates for i < 100 and contributes to 
50% or more of the total statistical error up to ^ ^ 200. 
Systematic errors due to uncertainties on the filter (in 
blue) and beam smoothing function (in yellow), which 
were computed as discussed in Sect.|21 are well below the 
statistical errors. 

Figure [7| shows the absolute value of the normalised 
error covariance matrix of the Archeops angular power 
spectrum for the binning shown in red on Fig.|Sl The cor- 
respondence between bin number and multipole range is 
indicated in Tab. ^ This matrix was computed using the 
simulations described in section rOl and provides the abso- 
lute correlation between multipole bins. The off-diagonal 



Fig. 6. Detailed description of the statistical error bars 
(in black) on the Archeops angular power spectrum ob- 
tained with Xspect in terms of sample variance (in cyan) 
and instrumental noise (in red). In addition, systematic 
errors on the angular power spectrum estimation due to 
uncertainties on the filter (in blue) and beam smoothing 
function (in yellow) are shown (see Sect. IS)). 




Fig. 7. Error covariance matrix of the Archeops angular 
power spectrum computed using the Xspect method. The 
correspondence between bin number and multipole range 
is indicated in Tab.^ The off-diagonal terms are less than 
12 % 
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Fig. 8. Archeops temperature angular power spectrum 
obtained using the SMICA method for one (in red) and 
two (in blue) components. A mixing of log- linear scales 
is presented to improve the readibility of the figure both 
on the Sachs-Wolfe plateau and on the acoustic peaks 
regions. 



terms are less than 12 %, and therefore the estimates 
can be considered as roughly uncorrelated across bins on 
multipole space. 

4.2. Archeops temperature angular power spectrum 
using SMICA 

To apply the SMICA method to the Archeops data we 
choose to estimate two components (number required by 
the data: see figure El and related comments) correspond- 
ing to the CMB anisotropics and to unidentified residuals 
from foregrounds. The mixing matrix is simultaneously es- 
timated allowing for recalibration of individual detectors 
against the most sensitive photometric pixel at 143 GHz. 

We find that CMB anisotropics are clearly detected for 
all the bolometers. A second component, much weaker in 
amplitude, is significant only in the 217 GHz maps. This 
component is thought to be a weak residual of foreground 
subtraction (see Sect. 15.31 for a more detailed discussion). 
FigurelSlshows in red the estimated CMB power spectrum 
with SMICA assuming two components. 

To assess the impact of the second component, we 
run SMICA assuming a single physical component in the 
Archeops maps, meant to be the CMB anisotropics. For 
this second analysis, we fix the CMB mixing parameters 
to the values derived from the dipole calibration, allowing 
the direct comparison with Xspect. Figure |S1 shows in blue 
the CMB power spectrum obtained in this way. 




800 



Fig. 9. Rescaled spectral fit as function of the multipole. 
The dashed line is the fit for 1 component, the solid line 
curve is for 2 components. The dotted lines are the bounds 
of the 68% confidence interval estimated in simulations 
of the two-component model. Note how a two-component 
model brings the spectral mismatch within the statisti- 
cal error bounds, showing that in addition to the CMB 
anisotropics a second component is required by the data 
mainly at low spatial frequency {£ < 100). 



The fit of the estimated model to the data is quantified 
by the lowest possible value (l){6) = mine 0(^') of the spec- 
tral matching criterion Eq. Q). If the model of observa- 
tions is correct (i.e. includes the probability distribution 
of the data), then <j){6) should be statistically small. A 
finer picture is obtained by splitting the overall fit of (j){9) 
into its components 'WqK{Rq, Rq{9)) as a function of the 
multipole bin q. Figure [HI shows the spectral adjustment 
of the best one-component model and of the best two- 
component model. The adjustment is much better with 
two components than with a single component, indicating 
that a second component is required by the data. 

Blind estimate for two components allows to separate 
systematic residuals in the two 217 GHz maps at the cost 
of some small increase in the CMB power spectrum error 
bars. The errors on the estimated CMB mixing parameters 
(bolometer intercalibration error) influence the error bars 
on the power spectrum estimate. The ratio between CMB 
power spectrum statistical error bars for the two and one 
component cases is about 20 % at low £ and 10 % at high 
i. 



5. Discussion 

The CMB angular power spectrum measured by Archeops 
as computed using Xspect and SMICA extends to 
a larger multipole range the results presented in 
IIBenoit et al. 2003a|l and is in good agreement with them 
on the common multipole range reducing the error bars 
by a factor of three. 



10 



M. Tristram, G. Patanchon, J.F. Maci'as-Perez et al.: Archeops CMB power spectrum 



6000 




-4000 



100 200 300 400 500 600 700 

I 



6000 [ 



4000 



+ SMICA 2c 
+ Xspect 



^ 2000 



-2000 ■ 



-4000 L 



9 f 




10 



100 200 300 400 500 600 700 



Fig. 10. Xspect angular power spectrum using six detec- 
tors (in red) compared to the one obtained using only the 
four 143 GHz detectors (in blue). The difference between 
the two power spectra is given in the bottom plot (shifted 
by 2000) and are compared to the error bars (black dotted 
line) . 



5.1. Consistency checks 

Internal tests of consistency have been implemented in 
order to check the robustness of the results presented 
above. The Archeops CMB angular power spectrum has 
been computed for two different map resolutions {nside — 
512,256 corresponding to 7 and 14 arcmin pixels resp.) 
and we observe no significant differences between them. 
Furthermore, we have substantially varied the frequency 
intervals for the timeline bandpass filtering and no signifi- 
cant effect appears in the estimation of the angular power 
spectrum even at high multipoles. In addition, to check the 
consistency of the results between the two CMB channels 
(143 and 217 GHz) we have computed, using Xspect, the 
CMB angular power spectrum for only the four 143 GHz 
bolometers. Figure 11(11 shows this spectrum (in blue) com- 
pared to the one using the 6 most sensitive photometric 
pixels (in red). The spectra are in very good agreement, 
within the error bars, over the full multipole range. Using 
only the 143 GHz bolometers reduces significantly the sen- 
sitivity to the second acoustic peak but no systematic off- 
sets are observed. 

As an extra consistency check, we compare in Fig. 1111 
the Archeops angular power spectrum obtained using 
Xspect (in red) with the one computed with 2-components 
SMICA method (in blue). The difference between the two 
power spectra, given in the bottom plot, is well below the 
error bars (red and blue dotted line). Detailed discussion 
of this issue is presented in Sect. 15.31 



Fig. 11. Archeops angular power spectrum using 
Xspect (in red) and using SMICA (in blue) . The difference 
between the two power spectra is given in the bottom plot 
(shifted by 2000) and are compared to the error bars (blue 
and red dotted line). See text for details. 
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Fig. 12. Fourier spectra obtained through the use of the 
two methods described in Sect. l3.5l for the best Archeops 
bolometer at 143GHz. These analyses are in agreement 
within the error bars. 



5.2. Archeops temperature angular power spectrum 
on the rings 

We show on Fig. El the Fourier spectra obtained through 
the use of the two ring analysis methods described in 
Sect. 13.51 for the best Archeops bolometer at 143GHz. 
These analyses are in agreement within the error bars and 
show a clear detection of the first acoustic peak. These 
results indicate that the processed timelines contain no 
obvious spurious feature at a particular time frequency. 
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5.3. Contamination from foregrounds 

As any balloon-borne experiment, Archeops is exposed 
to the fluctuations of the atmospheric emission. Moreover 
the Galactic emission at 143 and particularly at 217 GHz 
is low but not negligible. Even if a careful decorrelation 
to suppress ozone and dust spurious emissions has been 
performed (see Sect. I2.2|l . the residuals from this decor- 
relation are a potential source of systematic errors in the 
determination of the CMB angular power spectrum. 

The Galactic dust contribution must be much weaker 
at high Galactic latitudes. To assess the level of Galactic 
residuals, we have computed the angular power spectrum 
of the Archeops data using only the Northern part of 
the Archeops sky coverage. Figure El shows the es- 
timate of the angular power spectrum for the Galactic 
mask described in Sect. 12.21 (in red) and for high positive 
Galactic latitudes: b > -f20deg. (in blue) using Xspect. 
The differences between the two power spectrum esti- 
mates, shown in the bottom plot, are significantly smaller 
than the error bars associated to them. We conclude from 
this that the residual dust emission in the CMB angu- 
lar power spectrum obtained from the Archeops data is 
small compare to the statistical errors in the multipole 
range 17 < £ < 700. The multipole bin 10 < £ < 17 
shows a more important contimation from dust residual 
emission but still at the levels of the statistical and sys- 
tematic uncertainties. For £ < 10 we found that the dust 
contamination was significant and therefore this multipole 
range was not included in this paper. The same test has 
been performed using SMICA and leads to identical con- 
clusions. 

To fully assess the residual contamination to the 
Archeops data from dust and atmospheric emissions we 
have performed two independent tests based on Xspect 
and SMICA respectively. 

First, using the Xspect method we can cross-correlate 
the sky maps at 143 and 217 GHz used for the Ce estima- 
tion with the sky maps of the 353 GHz Archeops detec- 
tors. The observed emission on the latter is dominated by 
dust and atmospheric emission and to first order we can 
neglect the CMB emission. Thus from this cross correla- 
tion, we can obtain an estimate of the residual foreground 
contribution to the Archeops CMB angular power spec- 
trum computed with Xspect. The results from this anal- 
ysis are shown on Fig. 1141 The estimated contamination 
(in red) remains significantly below the statistical errors 
(in black) over the full multipole range except for the first 
multipole bin (£=[10-17[) for which the contamination is 
still smaller than the statistical error bar. 

As discussed in Sect. 14.21 we have performed, using 
SMICA, a two component analysis of the Archeops six 
best photometric pixels. The first component on this anal- 
ysis was identified as CMB emission whereas the second as 
the spurious residual foreground emission. This is signifi- 
cant only for the 217 GHz bolometer maps. This compo- 
nent is mainly due to residual atmospheric emission left 
behind after the linear decorrelation. This estimation is 
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Fig. 13. Xspect Archeops power spectrum computed 
for the Galactic mask described in Sect. 12. 2| fin red) and for 
6 > 20 (in blue) . The difference between the two estimates 
is given in the bottom plot (shifted by 2000) and compared 
to the error bars (blue and red dotted fine). 



represented in Fig. 1141 in blue, and can be compared to the 
foreground residual contamination estimated with Xspect 
at high multipoles. The SMICA estimate is of the same of 
order of magnitude and oscillates for £ > 200. These oscil- 
lations come from the uncertainties on the estimation of 
the second component which are well reflected on the error 
bars obtained for it. This could be due to correlated noise 
between the 217 GHz maps which would not be present in 
the residual foreground estimate obtained using Xspect. 
Further, this conclusion is reinforced by the fact that this 
contribution does not seem to be fully additive as expected 
from the SMICA model. 

From the above results we can conclude that the 
Archeops CMB angular power spectra obtained using 
Xspect and SMICA are fully compatible if we take into 
account the residual atmospheric contamination which is 
in any case well below the statistical error bars as shown 
in Fig. El We have plotted the differences between the 
Archeops CMB angular power spectra computed with 
SMICA 1 and 2 component (in blue), SMICA 1 compo- 
nent and Xspect (in black) , and SMICA 2 components and 
Xspect (in red). For comparison the statistical error bars 
are shown (black dashed line). This figure visually con- 
firms the fact that the contamination from foregrounds 
on the Archeops CMB angular power spectrum is well 
below the error bars. This analysis of the foreground con- 
tamination validates our choice of the galactic mask de- 
scribed in Sect. 12.21 
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Fig. 14. Atmospheric and Galactic dust residual emis- 
sions on the Archeops 143 and 217 GHz maps. In 
red, the residual foreground emission computed by cross- 
correlating these maps with the 353 GHz maps using the 
Xspect method. In blue, the residual foreground emis- 
sion obtained from the second component detected by the 
SMICA 2 components analysis of the Archeops data. 
In black, we plot for comparison the error bars of the 
Archeops CMB angular power spectrum. 



Finally, the contribution from point sources is negli- 
gible in the multipole range considered here (see paper 
I)- 

5.4. Comparison to the standard A-CDM cosmological 
model 

To check the validity of our results and their agreement 
with previous cosmological observations we have com- 
pared the CMB angular power spectrum measured by 
Archeops to the best-fit A-CDM cosmological model pre- 
sented in fSper gel et al. 2003| ). This model was derived 
from a combination of the WMAP data with other finer 
scale CMB experiments, ACBAR and CBI and is de- 
fined hy h = 0.71992, flth^ = 0.02238, flmh^ = 0.11061, 
T = 0.11026, constant n^ (0.05 Mpc)-^ = 0.95820 and 
normalization amplitude A (0.05 Mpc) = 0.73935. 

In Fig. El we present the best-fit A-CDM cosmological 
model described above superimposed on the Archeops 
CMB angular power spectrum which is rescaled by a factor 
1.07 in temperature (1.14 in Ci). This factor has been 
computed by assuming that the differences between the 
Archeops data and the model are due to a global scaling 
factor for all angular scales which has been fitted to 1.07± 
0.02 with of 27/24 and probabihty Q = 0.72. For this 



Fig. 15. Differences between the Archeops CMB angu- 
lar power spectra computed with SMICA 1 and 2 com- 
ponent (in blue), SMICA 1 component and Xspect (in 
black), and SMICA 2 components and Xspect (in red). 
For comparison the statistical error bars are shown (black 
dashed line). 



fit we have only considered the statistical error bars on 
the angular power spectrum. 

We observe that the agreement between the rescaled 
Archeops data and the model is very good. Here the 
model can be thought of as a guideline summarising other 
CMB experiments at different frequencies, in order to 
show the overall consistency across the electromagnetic 
spectrum. The scaling factor can be explained by the un- 
certainties on the absolute calibration of the Archeops 
data which are 6% in temperature (12% in Ci). A more 
detailed analysis of this issue is reported to a forthcoming 
paper including the determination of cosmological param- 
eters from the Archeops data as well as a comparison to 
other CMB observations at the map level. 

6. Conclusion 

Archeops was designed as a test-bench for Planck- 
HFI"^ in terms of detectors, electronics, cryogenics and 
data processing. Archeops has demonstrated the validity 
of these technical choices two years ago by determining, for 
the first time and in a single balloon flight, the tempera- 
ture angular power spectrum of the CMB from the Sachs- 
Wolfe plateau to the first acoustic peak (15 < £ < 350) 
using only two detectors. 

In this paper we present an improved analysis of the 
Archeops data using the six most sensitive detectors and 
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Fig. 16. The Archeops temperature angular power 
spectrum rescaled by a factor 1.07 in temperature super- 
imposed on the A-CDM best-fit model by the WMAP 
team and presented in ( |Spergel et al. 20031 ). 

20 % of the sky, mostly clear of foregrounds. Archeops 
has measured the CMB angular power spectrum in the 
multipole range from £ = 10 to ^ = 700 with 25 bins, 
confirming strong evidence for a plateau at large angular 
scales followed by two acoustic peaks centered around i — 
220 and I = 550 respectively. 

The Archeops CMB angular power spectrum has 
been determined using two different statistical methods, 
Xspect and SMICA. The results from these two methods 
are in very good agreement with differences between them 
well below the statistical error bars. Furthermore, they al- 
low a detailed analysis of the residual foreground contri- 
bution which is mainly due to atmospheric and Galactic 
dust emissions. The residual foreground emission on the 
Archeops data is small with respect to the error bars at 
all multipoles. 

Finally, we have compared the Archeops CMB angu- 
lar power spectrum to the best-fit A-CDM cosmological 
model presented in ( |Spergel et al. 2003| ) derived from a 
combination of the WMAP data with other smaller scale 
CMB experiments (ACBAR and CBI). We find that the 
Archeops data are in very good agreement with this 
model considering a rescaling factor to account for un- 
certainties on the absolute calibration. 

A more detailed analysis for the determination of cos- 
mological parameter with Archeops and other cosmo- 
logical datasets will be discussed in a forthcoming paper. 
Furthermore, a comparison of the maps from Archeops, 
WMAP and other CMB datasets will be used to study 
the primordial nature of the measured CMB anisotropics 
from their electromagnetic spectrum. 



All methods developed for this analysis will be im- 
plemented for the Planck-HFI data analysis. Even if 
Planck is less prone to systematic effects due to its space 
environment, the know-how acquired on Archeops data 
should prove useful in order to assess Planck final power 
spectrum. 

Acknowledgements. We would like to pay tribute to the mem- 
ory of Pierre Faucon who led the CNES team on this successful 
flight. The HEALPix package was used throughout the data 
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bin 


P ■ 


i-max 




XSPECT 

"totajl crroi' 


iiistj ruiTicnta,! error 




SMICA 

total error 


instrumental error 


sample variance 


1 


10 


16 


774 


251 


45 


899 


217 


11 


206 


2 


17 


24 


998 


167 


12 


1027 


170 


15 


155 


3 


25 


34 


1043 


168 


41 


999 


149 


22 


127 


4 


35 


49 


1487 


144 


39 


I486 


131 


26 


105 


5 


50 


59 


1217 


185 


51 


1081 


173 


39 


134 


5 


60 


69 


1537 


195 


54 


1561 


189 


48 


141 


7 


70 


79 


1613 


227 


78 


1619 


206 


57 


149 


8 


80 


89 


2038 


234 


78 


1978 


223 


67 


156 


9 


90 


99 


2275 


258 


93 


2451 


242 


77 


165 


10 


100 


119 


2586 


204 


74 


2639 


201 


71 


130 


11 


120 


139 


3193 


238 


90 


3221 


232 


84 


148 


12 


140 


159 


3148 


273 


110 


3234 


274 


111 


163 


13 


160 


179 


4225 


312 


138 


4358 


312 


138 


174 


14 


180 


199 


4941 


339 


159 


5050 


356 


176 


180 


15 


200 


219 


4589 


369 


189 


4506 


377 


197 


180 


16 


220 


239 


5085 


392 


219 


5183 


388 


215 


173 


17 


240 


259 


4258 


421 


263 


4340 


402 


244 


158 


18 


260 


279 


4356 


374 


235 


4538 
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226 


139 


19 


280 


309 


3174 


325 


233 


3385 


302 


210 


92 


20 


310 


349 


2325 


302 


247 


2351 


298 


243 


55 


21 


350 


399 


1960 


322 


292 


1862 


309 


279 


30 


22 


400 


449 


1832 


418 


394 


1825 


399 


375 


24 


23 


450 


524 


2569 


507 


483 


2487 


465 


441 


24 


24 


525 


599 


2394 


799 


774 


2649 


676 


651 


25 


25 


600 


699 


1885 


1183 


1168 


1595 


1124 


1109 
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Table 1. Archeops CMB power spectrum and statistical error bars (total, instrumental and sample variance) in 
(/iKcMe)^ computed with Xspect and SMICA (with two components) for the best six photometric pixels. 
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